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Abstract 

We show non-existence of solutions of the Cauchy problem in M. N for the nonlinear 
parabolic equation involving fractional diffusion d t u + (-A) s = 0, with 0 < s < 1 

and very singular nonlinearities <t> . More precisely, we prove that when (j>(u) = — 1 /u n 
with n > 0, or <j>(u) = log u, and we take nonnegative L 1 initial data, there is no 
(nonnegative) solution of the problem in any dimension N > 2. We find the range of 
non-existence when N = 1 in terms of s and n. The range of exponents that we find 
for non-existence both for parabolic and elliptic equations are optimal. Non-existence 
is then proved for more general nonlinearities <f>, and it is also extended to the related 
elliptic problem of nonlinear nonlocal type: u + (-A) s <f>(u) = / with the same type of 
nonlinearity <j>. 
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1 Introduction 


In this paper we discuss the question of non-existence of solutions of a class of nonlinear 
parabolic equations involving fractional diffusion and singular nonlinearities of the so-called 
ultra-fast diffusion type. These equations have the form 

(1.1) d t u + (-A) s 0(it) = 0, 

where 0 : M_|_ —> R is a monotone nondecreasing function of u with a singularity in u = 
0, namely such that 0(0+) = — oo. Consequently, nonnegative data and solutions are 
considered. On the other hand, (-A) s is the fractional Laplacian operator with 0 < s < 1, 
that is the nonlocal operator defined, at least for functions in the Schwartz class S, by 

(1.2) (-A) s v(x) = c(N,s)p.v. f Vx£R n , 

Jr* \x - y | JV +^ 

see, e.g., m- Here, p.v. stands for principal value of the integral and c(N, s ) is a positive 
scaling constant, whose exact value is not important throughout this paper. 

We consider the Cauchy problem posed in the whole space with initial data 

(1.3) u(x, 0) = uq(x) for x € M. N , N> 1, 

and we assume that uq is nonnegative and integrable, which corresponds to the natural con¬ 
ditions assumed for a mass distribution that evolves in time by nonlinear nonlocal diffusion. 
Accordingly, we look for nonnegative solutions that are integrable in space for every time 
t > 0. 

The limit case s = 1 of equation (|1.1D corresponds to the well-known model of nonlin¬ 
ear diffusion driven by the standard Laplacian, dtu = A 4>(u), sometimes called Filtration 
Equation. This equation has been studied by numerous authors in the last century, see the 
early works EDI and the monographs EZUGEj and their references. The well developed 
theory implies that given no as above, if 0 is continuous, strictly increasing and 0(0) = 00 
then there exists a unique solution of the Cauchy problem, where solution is understood in 
some suitable generalized sense like mild solution. Moreover, the operator that assigns to 
any initial condition its evolution is a semigroup of contractions in L^M^) (or in 
if two-signed solutions are considered). The set of admissible 0 includes, in particular, all 
powers 4>(u) = u m for m > 0, in which case the theory is more detailed. The use of mild 
solutions is usually changed into other familiar concepts, like weak or strong solutions, with 
better functional properties. Uniqueness theorems are proved to make such concepts useful. 

The extension of this theory to the fractional Laplacian model (11.1(1 has been done recently 
in a series of papers usi ns na eu for convenient choices of the nonlinearity 0. The first 
two references treat the case where <ft(u) is a positive power and construct a contraction 
semigroup with good properties, the third treats a case of logarithmic nonlinearity, while 
the last one treats the case of a more general smooth function 0, and then proves that when 
4>'(u ) > 0 for all u (i. .e., 0 is non-degenerate) the nonnegative solutions are indeed positive 

1 Actually, less stringent conditions on 0 are acceptable. 
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and C°° smooth. The quantitative analysis of positive solutions was then pursued in [5] 
in the form of regularity estimates, together with the study of existence and uniqueness of 
initial traces. 

Non-existence in singular diffusion. The purpose of this paper is to show a phe¬ 
nomenon of non-existence of solutions of the Cauchy Problem (ll.ll) - (ll.3f) when the non¬ 
linearity 0 is of singular type, more precisely when 0(u) is defined and increasing (and 
possibly smooth) for u > 0, but 0(0+) = — oo. This forces us to work only with nonneg¬ 
ative solutions, a restriction that is kept throughout in the paper. This type of singular 
nonlinearities were first considered for the Filtration Equation and then in a number of 
physical applications. The prototypical choice for the nonlinearity is an inverse power, that 
is 4>(u) = cu m for m = —n < 0. As an example, when m = — 1, the equation rules the 
evolution of the temperature in the Penrose-Fife model for phase transition (see |19] and 
pT] and the references therein). 

The first caveat is that one has to change the sign of the power, c = —a < 0 if we want 
the resulting singular equation to be formally parabolic, for then we can write dtu = A 4>{u) 
in the form 

(1.4) dtu = div (H(u)Vu), with D(u ) = > 0. 

Since the diffusivity D(u) is infinite at u = 0 we talk about singular diffusion or fast 
diffusion, and for n > 0 (i.e., m < 0) we call it ultra-fast diffusion or very singular diffusion, 
cf. [28j. There is no loss of generality in assuming that an = 1 to simplify the form of 
the equation. The second caveat is that the case m = 0 enters into this theory if we put 
0(u) = log(u), since then D(u) = l/u, cf. (28] and its references. 

Very singular fast diffusion equations have curious properties, like the non-existence of 
solutions for reasonable classes of data that we are going to prove. The present work is 
motivated by the result of [26], where it was found in that no solutions exist for the Cauchy 
Problem for equation (11.41) with nonnegative, integrable initial data in the following cases: 
if N > 3 for m < 0, or if N = 2 for m < 0, or if N = 1 for m < —1, and besides the 
range is optimal (note that we write m = — n). And the nonexistence applies to any time 
interval, however small, we have a phenomenon of instantaneous extinction or extinction in 
zero time. It must be noted that the problem with existence occurs exclusively because of 
the singular level u = 0. Thus, solutions of this problem do exist for data that do not go 
to zero at infinity, or even for data in some L P (M N ), on the condition of not having a rapid 
decay as \x\ —> oo. Another situation in which solutions exist is the case of the equation 
posed in a bounded domain coupled with zero Neumann conditions or even with dynamic 
boundary conditions (see [22]). 

One then wonders if the presence of fractional diffusion will affect the non-existence result, 
since in particular we cannot identify a diffusivity as the D(u) in the standard Laplacian 
case; if the answer is yes, we want to identify the techniques that enable us to treat the 
problem and get optimal results. We will give answers to such questions below. Briefly 
stated, the phenomenon of non-existence of solutions still holds, and we find that the best 
approach to the proof is via limits of non-degenerate regularizations and the use of Riesz 
potentials. Moreover, when 0(u) has the form of an inverse power with exponent n, we get 
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the optimal range of exponents n for every dimension N and fractional exponent s. And 
we extend the result to more general nonlinearities, though in a less sharp formulation. We 
also obtain sharp non-existence results for elliptic problems. 

The next section contains a statement of the problem, the technique of approximation by 
regularization, the concept of limit solution, and the statement of the main results with 
optimal ranges of application. Further outline of the contents of this paper will be given at 
the end of that section. 


2 Problem, limit solutions and main results 


We discuss the question of non-existence for the following class of diffusion equations 


d t u + (-A) s cf)(u) = 0, in Q, 
u(x, 0) = uq(x) for x € , 


where Q := R N x (0,T), with T > 0 and N > 1. The nonlinear function is defined, 
increasing and smooth for u > 0 with cj)(s ) —> —oo as s —> 0. More precisely, in this section 
and in most of the paper cj) chosen from the list 


( 2 . 2 ) 


4>n(u) 


—u n , for n > 0, 
log a for n = 0. 


We remind the reader that we will use throughout the notation N for the space dimension, 
and n for the nonlinearity exponent (m = —n is the more standard notation used when 4> 
is not singular, im but using m < 0 is rather inconvenient). The reason for which the case 
m = 0 is rewritten in terms of the logarithm nonlinearity is well-documented for classical 
diffusion s = 1, where the equation 


d t u = A(logu), 


can be rigorously seen as the limit as m \ 0 of the porous medium equation 

dtu + (-A ) s (u m /m) = 0, m > 0. 


cf. [20], [28( Chapter 8] and also [3]. For 0 < s < 1 similar results are obtained in [15] and 
in the recent paper [25] for N = 1. 

Following |26j , the strategy of proof of our nonexistence results is based on approximating 
problem (12.11) with data uq € T 1 (M iV ) by the family problems 


d t u £ + (-A ) s (4>(u e )) = 0, n> 0, 
u £ (x, 0) = uq(x) + e for x € R N , 


so that we avoid data with values on the singular level u = 0, so that any £ > 0 consider 
the sequence uo,e '■= uq + e. Then, the standard theory must apply, see next section, and a 
classical solution u £ (t,x) exists for all e > 0 and u £ > e. Moreover, the maximum principle 
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holds for these classical solutions and we have u £ > u £ > for e > e' > 0. Therefore, we will 
be allowed to take the monotone limit 

(2.4) u(x, t) = lim u £ (x, t ), 

£->0 

see the proof of Proposition 13. II below for a rigorous justification. This function is a kind 
of generalized solution of the problem, that we will call the limit solution. It is now an 
important step of the theory to decide in which sense this is a solution of the equation in a 
more traditional functional sense (like weak, strong or viscosity solution) and also in which 
sense it takes the initial data. In cases of non-uniqueness of such solutions, the unique limit 
obtained by the above method has been called by various names: maximal solution, SOLA, 
proper solution,... 

There is an extreme possibility that we will discuss in this paper: in the limit e —>• 0 the 
family u £ might converge to zero for any t > 0 and all x £ M. N , namely u(x, t ) = 0, and this 
is the plain statement of the non-existence theorem. For clarity of presentation, we split 
the results in the following two Theorems. 

Theorem 2.1. Let N > 2, n > 0, and 0 < s < 1. Then the limit solution u of m 
corresponding to any initial data any uq £ L l + {H N ) vanishes identically for all x £ , 

t > 0. More precisely, we have 

(2.5) lim u £ (x, t) = 0, 

£->-0 

and the convergence is uniform in W N x (r, oo) if uq is bounded, while it happens locally in 
L 1 (Q) for general uq £ L^R^). More precisely, for every ball B C R^ 

(2.6) lim u E = 0 in L 1 {B) uniformly in t > t for any r > 0 . 

£\0 


The one-dimensional case is worth a separate statement since there appears a critical 
relation, 2s — n = 1, to determine the range of existence or non-existence. It will require a 
special analysis. 

Theorem 2.2. Let N > 1, n > 0, and 0 < s < 1. Then the limit solution u of m 
corresponding to any initial data for any uq £ -^(R) vanishes identically in the sense of 
the previous Theorem unless 1/2 < s < 1 and n < 2s — 1, or s = 1/2 and n = 0. This range 
of exponents is optimal. 

The non-existence picture in one dimension is given in Figured] 

Remarks, (i) It is important to notice that the function u = 0 in R^ x (r, +oo), Vr > 0 
cannot be a solution of the equation in any differential sense since (f(u) = —oo for any point 
of R^ x (r, +oo), Vt > 0. On the other hand, the initial value has collapsed from uq{x) + e 
to zero at t = 0+, and not to uq(x) as expected. This is called instantaneous extinction, 
and shows a strong discontinuity of the semigroup at t = 0, for data rio £ L^R^). 

(ii) The result depends in a crucial way on the class of initial data. Notice that positive 
constants are indeed classical solutions to this problem and of course they do not extinguish. 
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Figure 1: The non-existence picture in one dimension. The white upper triangular 
region, corresponds to the range of parameters where there is existence: n < 2s — 1 and 
1/2 < s < 1, together with the special point n = 0 and s = 1/2. The grey region corresponds 
to the range of parameters where there is no existence. 

Therefore, there should be an intermediate decay between L 1 data and constants, such that 
nontrivial solutions with that decay at infinity will still exist. For standard diffusion s = 1, 
the critical decay that implies existence has been investigated by Daskalopoulos and del 
Pino in the series of papers mmm in which they provide an optimal decay condition 
to determine the existence or the immediate extinction in classes of L p data. 

(iii) Technically, the proofs of the above theorems require a different treatment, according 
to whether N > 2s or not. The critical situation happens when N = 2.s, hence there appear 
two different types of behaviour only in the one dimensional case. Indeed, for all dimensions 
N > 2 we always ve JV >2 s and this is covered in Proposition 14.11 and the same proof 
allows also to show non-existence the case N = 1 and 0 < s < 1/2. 

Once we have proved that the approximate solutions converge to zero, the natural question 
is what happens to solutions in a more traditional sense, like weak, strong or viscosity. To 
answer this question, we first we have to introduce the concept of solutions that we are 
going to use. 

Definition 2.1. A function u is a weak (energy) solution of equation (11.11) if it is nonnegative 
and 


• u € C°(0, + 00 ; L 1 <j>(u) <E Lf oc (0, + 00 ; H s (M n )) 


• it satisfies, for all C £ Cq°(Q), 

POO P POO P 

(2.7) / / udtCdxdt — / (-A) s / 2 <f>{u) (-A ) s ^ 2 (dxdt = 0. 

Jo Jr n Jo Jr n 

• A weak energy solution is said to be a strong solution if moreover u £ ^((O, + 00 ); L 1 (M iV )) 
and dtu(t, •) € L 1 ( R^) for every t > t > 0 and every r > 0. 


5 






We next establish a comparison result between approximate and strong solutions, namely 
we show in Theorem rm that u < u £ where u is a strong solution and u £ is the above 
mentioned approximate solution. More details are given in Section [7J where we prove the 
following theorem. 

Theorem 2.3. Let u denote a nonnegative strong solution of m, in the range of pa¬ 
rameters ( N,s,n) allowed in Theorems 12.11 and \2.2\ Then, u = 0 in M. N x (r,+oo) for any 
t > 0. This range of exponents is optimal. 

We shall mention that once we establish our non-existence result for strong solutions, 
then the same result will hold for any reasonable limit of strong solutions. Moreover, the 
Theorem above entails that no Barenblatt solution exists for the singular equation (12.11) . 
More generally, we can conclude that there exist no solution (in a reasonable differential 
sense) to m that regularizes into L 1 (R N ) for positive times. 

Outline of the rest of the paper. Further results. We establish the existence of 
solutions of the approximation problems and their properties in Section [3l and pass to the 
limit to obtain a so-called limit solution. 

We then prove the non-existence results, in the form of vanishing limit solution, first for 
N > 2s, Section 0J and then for the remaining cases in ID, Section [5l We show optimality 
of the exponents in ID in Section [6j Sections 01 [5] and [6] contain the proof of Theorems 12.11 
and 12.21 

Section 7 contains the proof of above mentioned comparison result between strong and 
approximate solutions, together with the proof of Theorem 12.31 

In later sections we deal with two interesting extensions of the results. 

• More general nonlinearities. Using concentration comparison and rearrangement tech¬ 
niques developed in [301 EQ, we prove that the results of Theorems 12.11 and 12.21 can be 
extended to more general nonlinearities, as stated in the Theorem 18.21 we refer to Section 
[8] for more details. 

• Elliptic equations. An important extension of our results concerns the question of 
existence for solutions to elliptic equations of the form 

(2.8) u + (-A) s (f{u) = / , with / e Zq_( W N ). 

A complete study is performed in Section 0] where we show how the parabolic technique 
developed in Sections 01 [5] and [6] can be adapted to treat this problems and obtain analogous 
non-existence results, that are stated in Theorem 19.21 In particular, the optimal ranges of 
non-existence of the elliptic and parabolic problems coincide. 

There is another section devoted to the Dirichlet problem on bounded domains. This is 
rather straightforward, once we make sure that a suitable comparison holds. 
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3 Study of the approximate solutions 

It is convenient to write u £ (x,t) = v e (x,t) +s and then try to solve the Cauchy problem 


d t v + (-A) s (0 e (u)) = 0 with cj) £ (v) := 4>(v + e) - 0(e) 
n(0) = uq for x € M N . 


for all e > 0. This is a modified problem prepared to avoid the degenerate level of the 
equation by displacement of the axes. Note that for e > 0 and for nonnegative arguments 
c i> e is a smooth, positive, monotone increasing function with 0 e (0) = 0. 

Let us construct the approximate solution u £ to (12.31) . It will be convenient in a first 
approximation to assume that uq is bounded. 


Existence and properties. The theory of existence and uniqueness of weak solutions 
to (13.11) is given in [231 Theorem 8.2]. For any e > 0, let v £ denote the weak solution to 
m- Here weak solution means a function v £ C°([0, +oo); L 1 (R JV )) such that 0 £ (u) £ 
L'lod0,+oc; H s (M. n )) satisfying u(0) = uq a.e. in M. N and 

roo r roo r 

(3.2) / / vdtCdxdt- / (-A) s/2 (0 £ (u)) (-A) s/2 (dxdt = 0 

Jo Jr n Jo Jr n 

for any 0 £ C%°(Q). After obtaining these solutions we restore for any e > 0 the original 
u-level by defining u £ := v £ + e, as stated at the beginning. Clearly, we have that u £ > e 
(actually, u £ > e) and hence v £ > 0 in Q. 

Let us list some further properties of v £ and u £ . 

• Boundedness and regularity. These solutions are shown to be bounded for strictly positive 
times. More precisely, for every t > 0 and every p £ [1, oo] there holds 


(3-3) ll l, e('d)llLP(R iv ) ^ II^oIIlp^) 

As a consequence, if uq belongs to L°°, then v £ is regular enough to satisfy the equation in 
the classical sense at least when t > 0 , by the results of [24]. Therefore, u £ = v £ +e is smooth 
and satisfies the original equation in the classical sense in Q. Under these circumstances, 
the initial data are also taken, at least in the sense of convergence in L 1 (R JV ). 

• Conservation of the L°° bound. If u o belongs to L°°, then we can remove the boundary 
layer at t = 0. More precisely, we can proceed as in mm to conclude that any classical 
solution of (12.31) and of (13.11) satisfies the following 


(3.4) 


IM0IIl°°(iA) < vo < r < t. 


• Mass conservation. Nonnegative solutions to the evolution equation (13.11) conserves the 
mass, cf. mm~ More precisely, we have for all t > 0 


(3.5) [ v £ (x,t)dx= [ 

Jm. N Jr 1 


uo(x)dx 


i. e., 


/ 

Jr 1 


(u £ (x, t ) — e)dx 
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• L 1 -contraction and comparison. The evolution (13.11) is an L 1 contraction, namely we have 

(3.6) f ( u £ (t ) - z e (t))+ dx < [ (n 0 - z 0 )+dx for t > 0, 

Jr n J r n 

where is a solution to (12.31) constructed starting from another data zq G T 1 (M Ar ). Here 
and in what follows (■)_)_ denotes the positive part function. For the proof we may again 
refer to mm- 

• Monotonicity with respect to e. An easy version of the above comparison argument 
shows also that for 0 < e < e' we have 0 < e < u e < u e >. 

• Time monotonicity. The scaling invariance of the equation entails the following important 
mononoticity property of the strong solutions, known as the Aronson-Benilan inequality, 

u 

(3.7) d t u e < g . V(x,t) G Q. 

(n + l)t 

Proof. This is rather classical in the PMR theory, m , but we give here a proof for 
convenience of the reader. This inequality is a consequence of the scaling properties of the 
equation. We give a brief sketch of the argument that lead to (13.71) (for more informations 
see EH). Starting from a smooth solution u of (12.ip and a parameter A > 1, we construct the 
function u\(x, t ) := A u(x, A ~ n t). The scaling invariance of the equation gives that u is again 
a solution of (12.11) with initial condition u\(x, 0) = Arto(x). In particular, note that since A > 
1 , we have that u\(x, 0) > uq(x) and thus, thanks to the maximum principle, we have that 
u\ > u in Q. An important consequence of this fact is that the right derivative > 0 

(this inequality can be proved by working directly on the incremental quotient). Thus, we 
have 

0 < u\ = u — (n + 1 )tdtu. 

dA|A=i 

□ 

Passing to the limit. We may now pass to the limit e —> 0 using the monotonicity of the 
family u e . When uo is bounded the limit is taken in a set of uniformly bounded functions, 
hence by the monotone convergence theorem the limit u is taken in the local L 1 sense, 
i.e., in L l (B) for every compact subset of x [0, oo]. When uq is not bounded but still 
integrable, we use the L 1 (M Ar ) contractivity result to obtain the same result. Summing up, 

Proposition 3.1. If uq is a nonnegative function in L 1 (1R* /V ) there exists the monotone 

limit u = lim u £ with local convergence in L 1 ((5). 

£—>0 


4 Proofs of non-existence. The main case 

We will establish here the trivial limit of the approximation scheme in the range N > 2s > 0, 
of course for nonnegative data uq G L 1 (M' v ) and exponents n > 0. This subsection is the 
core of our technical argument. An important feature in our proof is that the inverse 
operator of (-A) s can be expressed as a convolution with the Riesz potential (see e.g.|23]), 



namely 


( 4 -l) (-A)-‘u(x)=c N ,,^ N ^-§-dy, 

where cn, s is a positive normalizing constant. The above representation does not hold in 
the subcritical case N — 2s < 0. This happens in dimension one and for s € [1/2,1), and 
requires a different argument, cf. Section [5l 

We want to show that the sequence u £ of solutions to (12.31) converges to zero as e \ 0 for 
strictly positive times. The following result sums up the situation. 

Proposition 4.1. For any given uq > 0 with uq £ T 1 (M Ar ), let u £ be the solution of (I2.3|) . 
Then for all n > 0 and N > 2s we have 

(4.2) lim u £ = 0 

and the convergence is uniform in M. N x (r, oo) ifuo is bounded (the time r > 0 is arbitrary), 
while it happens locally in T 1 (Q) for general uq £ L^R^). More precisely, for every ball 
B CR N 

(4.3) lim u £ = 0 in L l (B) uniformly in t > r for any r > 0 . 

e\0 

Notice that this result gives the proof of Theorem 12.11 i.e. applies to dimensions N > 2 
for all 0 < s < 1 and all n > 0. However, in IV = 1 we only cover the case 0 < s < 1/2, 
therefore we prove Theorem 12.21 only when 2s < 1. The whole picture in ID depends on 
other arguments and will be treated in the next section. 

4.1 Proof of Proposition 14. ll 

The proof is divided in several steps. We will begin by considering the case of initial data 
that are smooth and with compact support (or decay fast at infinity). Finally, in a the last 
step, using an approximation/density argument we will cover the general case of L+(M. N ) 
initial conditions. 

• Step 1. Proof of convergence for good data. Approximate solutions. Consider the 
sequence u £ of solutions to (12.31) . Recall that 

(4.4) e < u £ in 1^ x (0, +oo), 

(4.5) / (u £ (x,t) — e)dx = / uo(x)dx, Vi > 0, 

Jr n Jr n 

(4.6) the sequence u £ is monotone w.r.t. e 

From these properties we already know that have that the limit function u(x, t ) := liim^o u £ (x, t) 
is well defined and finite for any (x,t) £ Q. Moreover, we have that the function u is such 
that u(-,t ) £ L 1 (R Ar ) for any t > 0. More precisely, we have 

(4.7) l|H( i )llL 1 (K iv ) < II' u o||l 1 (r jv ) ^ > 0- 
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The main question is whether or not the mass conservation law is preserved in the limit, 
since both possibilities are compatible with Fatou’s theorem. We are going to prove that 
actually, mass is not conserved and the limit function is just zero. 

Let us take an arbitrary r > 0. We are going to prove that 
(4.8) u £ —°> 0 uniformly on M. N x (r, +oo). 


The uniform convergence is a consequence of the fact that the initial condition uq is assumed 
to be smooth and with compact support. We first need to prove the following lemma. 

Lemma 4.2. Under the running assumptions, the following equality holds true for any 
e > 0, 0 < t* < t and for any x £ 

(4.9) f <j){u £ {-,t))dt- {t-t*) 4>{£) = f U£<yV '^_ J^- V 2 S T W 

j t* j ~r n y i 


Proof. Recall that u E = v e + e, and that v e satisfies equation (13.11) . We will work with v £ , 
and we rewrite (14.9 j) in terms of v £ in the form: 


(4.10) 


[ (fe(v e (-,t))dt = [ 

Jt* Jr 1 


v e (y,T) -v e {y,r*) 


\x - y\ 


N-2s 


dy. 


The proof of this important identity can be obtained in several ways. One way is to consider 
in the weak formulation of equation (13.11) the following test function 


ips(t,x) 


X[r*,r]{t ) (-A) 


Xb 6 (o)( x ) 

1^5 (0)1 


Notice that for almost every x € M N , 


/ a \-s Xb s ( x)(y ) <5\o C 

\B s (x)\ \x-y\ N ~ 2s 


for almost every y £ W N 


Unfortunately, this function is not an admissible test function. One possibility is to approx¬ 
imate if>s by means of smooth functions, as it has been done for example in Proposition 4.2 
of [6j. Another possibility is to observe that v £ are smooth functions by construction, cf. 
[24 ] . in particular they are classical solutions to equation (13.11) therefore we can multiply 
the equation by ips obtaining 


(4.11) 


f [ ipsd t v e (y,t)dtdy = - [ f i/j s (-A) s 4> £ (v £ )dt dy 
Jr Jr n Jr Jr n 

Now we analyze the two terms separately. As for the left-hand side: 

Xb s ( o)(*) 


(4.12) 


[ [ i/j s d t v £ (y,t)dtdy = [ X[ T *,r](t) [ d t v £ (-,t) (-A) 

JrJr n Jr Jr n 

f r / \ / *',]/• X. (x) (jj) , (5N0 f 

= [vdy,T)-v e {y,T )](- A) dy - > 

Jrn \ b 5{x)\ J r i 


\B S m 

Ve(y,T) - v £ (y,r*) 

\x - y I 


dtdx 

N—2s d y ■ 
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the last limit can be justified by splitting the integral as follows. First of all, we notice that 


(4.13) 


(-A) 


- a XB S (x){y) 


< c- 


1 


x - y | 


N—2s 


m 


r>N 


for a constant C > 0 depending on N and on s. Then, we have that 
(4.14) 


Ve(y,r)-v E (y,r*) = f v e (x-y,T) -v e (x-y,T*) 

JR N /V 9 c 


/ 

J 


\x — y\ N 2s 
Ve(x — y,r) — v £ (x - y,r*) 


/ -Bi(O) 12/1 

< ||u e (-,r) -v £ (-,t*)\\ l . 


N-2s 


d y + 


\y\ N -2 S 

v e (x-y,T) -v e (x -y,r*) 


'flf(O) 


d y 


\y\ N 2S 

v £ (x-y,r) -v e {x -y,r*)\dy 


^i(o) \y\ N 2s y + Ibk o) 

< 2(C||u £ (0)|| Z/ cx 3 ( R jv) + ||t' £ (0)||ii(]KiV)) < 2(7 max{||uo||i,i(RiV) , ||«o||z,°°(rjv)} 


In the third line, we have used that is integrable in -Bi(O) and v £ (-,t) — v £ (-,t*) is 

bounded, while for |y| > 1 the kernel | y\ 2s ~ N > 1 and v £ (-,t) — v £ (-,t*) is in L 1 (R N ). We 
conclude then by dominated convergence. 

The right-hand side becomes: 


(4.15) 


i>5 


l 


(-A) s (f £ (v £ )dt dy = — f X[r*,T](t) [ (~ A ) S twttt Y^Y M v Y dt d d 

Jr 1 Jr* | 5 «(®)| 

J J>e{v £ {y,t))dt dy J (f> £ (v £ (x,t))dt 


I B s(x)\ Jb s (x) 

where in the second line we have used Fubini’s Theorem as follows: 
(4.16) 


Xb s {x)(z) ( f (-A Y 4> £ (v £ (y,t)) 


r xb S (x)(z) (r 
Jr n \ B 5(x)\ \Jr> 


\z - y\ 


[ (/ 

JR N \J K 

N- 2 s d y) dz = 

J Jr 


Xb s (x)( z ) dz 


r n \B s (x)\ \z — y\ N 2s 
Xb s (x) 


(-A) s (j) £ (v £ )dy 


\B S (x)\ 


(-A)~ s (-A) s Mve) d y 


In the last step we have used the definition of Lebesgue point for x i->- J^ (f> e (v £ (x,t))dt. 
Joining (14.121) and (14.161) we finally obtain (14.171) and hence (14.91) follows. □ 


• Step 2. Convergence as e \ 0 for bounded initial data. Now may now proceed with the 
proof that u £ ‘'^'°> 0, uniformly on R^ x (r, +oo). To this end, let us fix r and t* such 
that 0 < r* < r. The bound obtained in (14.141) combined with (14.91) . gives 


PT PT 

/ <fi(u £ {x,t))dt =— (f(u £ (x, t))dt > — (r — r*)<^(e) — C, Vx G 

Jt* Jt* 


t>N 
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where the constant C depends only on s and on N and is independent of e. As a consequence, 
if we let e \ 0 and recall that 0(e) ——oo (see (12.211 1. we have that 


(4.17) 


lim 

e\0 


f (j>(u e (x,t))dt 
Jr* 


= — oo, uniformly for x £ 


t,N 


Next, we use time monotonicity. Since, 


dtit-V^vit)) = t- l ^ n+1 \d t v{t) - —«(t)) Vt > 0, 

the validity of the Aroson-Benilan estimate (13.71) is equivalent to the decreasing monotonicity 
of t^ l ^ n+1 ^v{t) with respect to time. Thus, we have that u £ verifies 


(4.18) 


/ t \ XT 

u £ (x,t) < ( — ) u £ (x, s ), for any x £ and for any t > to > 0. 

\t 0 y 


On the other hand, using the definition of 0 in (12.21) for n > 0, monotonicity implies 


(4.19) 4>(u e (x,t)) < 


t \ -n/{n+l) 


to 


4>(u £ (x, to)), for any x £ M N and for any t > to > 0. 


Therefore, for any x £ M JV and n > 0 we have 

J <t>{u e {x, t))dt > C(t , t*, n)u £ (x, r)~ n , 

where the constant C can be computed explicitly. As for for n = 0, monotonicity implies 

(4.20) 0(ii £ (x,t)) < 4>(u £ (x, s)) + log , for any x £ R N and for any t > s > 0, 

Therefore, for any x £ and n = 0 we have 

[ (j>(u £ (x, t))dt > log(u £ (x, r)) - C(t, t*) 

Jt* 

where again the constant C can be explicitly computed. 


• Step 3. Uniform convergence for bounded data. All the above estimates imply that 

0(u £ (-,r)) —oo uniformly on M. N , which implies that u £ (-,t) —^ 0 uniformly in M w . 

Thus, using the L°°-decay (13.41) . we can conclude (14.81) . 

• Step 4. The case of unbounded initial conditions. When dealing with initial conditions 
which are only integrable, we will obtain a weaker convergence than (14.21) (as in the state¬ 
ment of Theorem 14.II) . by a standard density argument that we will discuss here. Given uq 
we prepare a sequence of bounded functions 0 < «o,M < M, and with compact support, 
such that 

uq,m ^ I — +00 > u 0 in 70(11^). 
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In particular, for any k > 0, we find M such that 

(4.21) \\uq,m ~ wo||li(r^) < p VM > M. 

Let u e> M be the solution to m starting from «o,m + £■ and let M > M. Since uo t M £ 
L°°(R n ), we know from the analysis above that in the appropriate range of N, n, s, we have 

lim u eM = 0, 

e\0 

uniformly in R^ X (r,+oo), Vr > 0. Fix now r > 0 and k > 0, and let M as in (14.211) . 
Then, there exists £o such that for e < e 0 , there holds 

UeAi(x,t) < — in R^ x (t,+ oo). 

’ k 

On the other hand, using (14.211) and (13.61) . we have for M > M 

I \Ue,M{-,t) ~ U £ (-, ^)HlI(rJV) < -. 

Thus, for every ball B in R^ we conclude that 

lim u £ = 0 in L l (B) uniformly in t > r for any r > 0 . 

e\0 

This concludes the proof of Proposition 14. 11 hence of Theorem 12.II and of Theorem 12.21 only 
when N > 2s. 


5 Trivial limit solution in dimension 1 

We now deal with the one-dimensional case, N = 1. This section contains the proof of 
Theorem 12.21 in the cases not included in the previous section (2s < 1), namely we consider 
here the remaining two cases: s € (1/2,1) with n > 2s — 1 and s = 1/2 with n > 0. 

Note also that these conditions are in complete agreement with the existing results available 
for the standard Laplacian, i.e. s = 1, where existence of nontrivial solutions holds when 
N = 1 and n < 1. Here we prove non-existence in the remaining cases. 

These two cases require a quite different technical treatments, mainly because the Green 
function has really a different form from the other cases, indeed it grows at infinity. 

The proof of the optimality of the range of exponents for this one dimensional case will 
be given in Section [6j 

5.1 Non existence for s £ (1/2,1) 

We now analyze the trivial limit in the situation N = 1, s £ (1/2,1) (so that the rule 
N > 2s does not apply) and n > 2s — 1. Thanks to the monotonicity (14.61) . the limit 

u := lim u E 

e\0 
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is well defined and finite for any (. x,t ) G Q = R x (0,+oo) and is such that u(-,t ) G L 1 (M) 
for any t > 0. We need a preliminary lemma. 

Lemma 5.1. Under the running assumptions, the following equality holds true for any 
£>0 ; 0 < t < f and for any 

(5.1) f <j)(u £ (x,t))dt= [ f(u £ (0,t))dt + [ p £ (-y)(\x+ y\ 2s ~ 1 - |y| 2s_1 )dy. 

J T J T J R 

where p £ (x) := u £ (x, t) — u £ (x, t), for i£l, 

Proof. For any 5 > 0, let Q denote a smooth function with compact support converging to 
the delta function in the origin as <5 \ 0, namely 

lim / fs(x)4>(x)dx = </>(0), smooth. 

s \°Jr 

We can assume that ^ > 0 and that f R (s(x)dx = 1. Then, denote with G$ the function 
given by expression 

(5.2) G s {x) := f \x- y\ 2s ~ 1 Cs{y)dy. 

J K 

Note that, even if 2s > 1, the function Gs is well defined and finite for any x in M due to 
the fact that (s has compact support. As a matter of fact, we have that (-A) S G <5 = C<5- 
Moreover, since f R (sdx = 1, using the elementary inequality 

(5.3) \A a -B a \<\A-B\ a , a € (0,1), VA, B G (0, +oo), 
we have 

(5.4) \G s (x) - M 2 *- 1 ! = | [ (Ix-yl 23-1 -\x\ 2s ~ 1 )Cs(y)dy\< [ \y\ 2s ~ 1 Cs(y)dy VrGl. 

J M J ]R 

Then, if we let S 0, we immediately get that Gs{x) ^°> |x| 2s_1 uniformly in R. Now, for 
x G R, we multiply the equation (12.31) with G$(x — •) — Gs(-)- and we integrate on Rx(r,f), 
with 0 <t <t. We get, 

(5.5) j p £ {y) (g 5 {x -y)~ Gs(yfj dy + ( W £ (y) - (t- r)^(e)) (( s (x -y)- (s(y))dy = 0, 

where 

(5.6) W £ {x) := J (j>(u £ (x,t))dt, for x G R. 

Note that since n > 2s — 1 > 0, we have <p{u £ ) = —uf n . Now, using (15.41) we have that, for 
any fixed x G R and for a sufficiently small 5, 

I Pe(y) (g$(x -y)~ Gsiy)^ \ < \p e (y)\ + \p e (y)\( I® - z/| 2s_1 - |y| 2s_1 ) I 

< (1 + \x\ 2s ~ 1 )\p e (y)\, Vt/GR 
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where we used (15.31) in the last inequality. On the other hand, we have that for any i£l 
the function y H > (1 + \x\ 2s ~ 1 )\p e (y)\ is in L 1 (M). Hence, recalling that (for fixed e) W £ is a 
smooth function, we can pass to the limit as 5 \ 0 in a standard way and obtain 

(5.7) W £ (x)-W £ ( 0)= [ p £ (—y)(jx + y| 2s_1 — \y\ 2s ~ 1 )dy , 

J i 

which is exactly m- □ 


We are now in the position to prove that u := lim^o^e = 0. We reason by contradiction 
and we assume that u( 0, t) > 0 for some t > 0. We keep the notations (15.61) for W £ and p £ . 

Combining (15.11) with (15.31) (in the equivalent form (15.71) . we obtain that W £ is Holder 
continuous, namely 

(5.8) | W £ (x) - W e ( 0)| < C\x\ 2s ~ l for iGl, 


with 

sup / |/O e |dx < C < +oo, 

£>0 Jw 

which gives 

W £ (x) > W e (0) — C'|a:| 2s - 1 . 

Thus, using the Benilan-Crandall estimate in the form (I4.19p . we have that there exists a 
constant C = C(r,i,n, ||'Uo||.l 1 (R)) such that 


(5.9) (j)(u £ (x,T)) > c(^)(u £ ( 0 ,t)) - \x\ 2s ^ for any 

Recaling that = —u~ n , we have that 

C 


x G 


1/n 


(5.10) 


U £ (x,T) > 


1 I I rp I 2s—1 

.u e ( 0 ,t) n ^ 


Thus, if we let e \ 0, the inequality above is conserved in the limit, namely 

(5- 11 ) u ( x , t )>(— T - ) • 

V s(TIF + l a; l 

Now, recall that we are assuming that u(0,t) > 0 (hence l/u(0,t) is hnite). This is man¬ 
ifestly impossible if n > 2s — 1 since the function (1 + |x|)~^ 2s_1 ^ n is integrable only if 
n < 2s — 1, while the function u(-,t) is in L 1 (M) for any t > 0. 

This concludes the proof of Theorem 12.21 when s G (1/2,1) with n > 2s — 1. 


15 








5.2 Non existence when s = 1/2 and n > 0 

Now we are going to discuss the case s = 1/2 and n > 0. Recall that for n = 0 there is 
an explicit positive solution, given in (16.21) . This situation presents some extra difficulties 
because for s = 1/2 the Green function is of logarithmic type and the corresponding integral 
analogous to the one in the right hand side of (15.71) seems difficult to handle. We rely on the 
analysis for s > ^ by noting that the lower bound (I5.10p . which is one of the main points in 
the argument above, has a “good dependence” on s, namely the constant C in (|5.10j) does 
not depend on s. Thus, our strategy in proving a lower bound like (15.101) for the solution of 
the s = 1/2 case consists in approximating this solution with solutions of the problem (12.31) 
with fractional laplacian of order s > 1/2. To make this argument rigorous, we need the 
following lemma on the continuity of the solution operator with respect to the fractional 
order of derivation 


Lemma 5.2. Let if : R —>• R be a smooth and non decreasing function with if(0) = 0. Let 

ss be a sequence of real numbers with s$ > ^ and such that s$ 1/2. Let v$ be the 

solution to 


(5.12) 


d t V8 + (-A) s,s tj’(vs) = 0 in Q 

u^(0) = vo t s in R. 


Then, if vo : s ^°> In T 1 (M), there holds that vs ^°> v in (^([O, +oo); L 1 (M)) and v is 

the solution to 


(5.13) 


dtv + (-A) 1//2 if(v) = 0 in Q 
u(0) = vo in R. 


Proof. The proof is an adaptation of the one given in [Til Theorem 10.1, Section 10]. □ 


We use the above Lemma in the following way. For any e > 0 denote with u e the solution 
to 

j d t u £ + (-A ) 1/2 {(j)(u e )) = 0, in R x (0, Too) 

\u e ( 0) = u 0 , e = u 0 + e in R. 

The above Lemma shows that u e can be approximated when 5 \ 0 with the solution us lS 
of 

f d t us >£ + (-A) S5 (4>(us, £ )) = 0, in R x (0, Too) 

1^(0) = u 0 ,s,e in R. 

To be precise, the above Lemma shows that the solution v e = u £ — e can be approximated 
by v s , £ := us, £ - £■ 

We are in the position to use Lemma 15.11 and repeat all the subsequent argument, for 
the approximations u§ £ which fall in that range of parameters. Denote, for any (x,t) £ 
R x (0, Too) , u(x,t ) := lim^o u £ (x,t) and assume by contradiction that u(0,t) > 0. Since 
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sg \j 1/2 as 6 —>• 0, we can fix 5 small enough in such a way that n > 2s^ — 1 > 0. Thus, 
we have that us >e satisfies the lower estimate (15.101) which becomes 


U8,e( x , r ) > 


1 . c 


1/n 


for some r <t. It is important to notice that the constant C that appears above, as well 
as the other constants that will come up in what follows, does not depend on s and hence 
on 5. Thus, if we let 6 \ 0, we have a similar lower bound for u e , namely 


/ C \ i/ n 

u e (x, t) > (-j-1 , for any i€l. 

V W + 1/ 


Now, we can let e 0, obtaining that 


f C \ 1 / n 
u(x, t ) > ( — y- 1 , for any i£l, 

' r.tn An + 1 ^ 


1 ( 0 , ty 


which contradicts the integrability of Thus, we have that u = 0 in R x (0 + oo). 

This concludes the proof of Theorem 12.21 when s = 1/2 and n > 0. 


6 Existence of nontrivial solutions in ID 

Let us now examine the situation for N = 1 in the remaining cases, namely 2s — 1 > n > 0, 
and s = 1/2, n = 0, where we expect to find nontrivial limit solutions, as announced in 
the statement of Theorem 12.21 This will prove the optimality of the range of exponents 
considered for non-existence and conclude the proof of Theorem 12.21 


6.1 Case s = 1/2, n = 0 


This is a kind of exceptional case in the parameter diagram, see Figure 1. Indeed, for 
s = 1/2 there exists an explicit solution of the evolution problem 


dtu + (-A) 1 ^ 2 (logu) = 0 inMx (0,+oo), 
u{ 0) = uq in M, 


with an initial condition in L 1 (M). It is given by the formula 

(6.2) U(x,t) = — —:—t^ in M x (0,+ 00 ). 

1 + \x\ z 

We see that it has the separate-variable type as in the previous subsection, and U(-,t ) is 
in L 1 (M) for any t > 0. It is very peculiar that the solution becomes identically zero in 
finite time. This is the so-called finite time extinction phenomenon which is typical of some 
ranges of fast diffusion, see ||28j for standard diffusion and mm for fractional diffusion. 
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A simple comparison theorem implies that any initial data uq(x) > c/(l + |x| 2 ) produces 
a nontrivial limit solution equal or larger than U(x,t ) with T < c/2, hence non trivial in 
a time interval 0 < t < T. In conclusion, by proving the existence of nontrivial solutions 
when n = 0, we have made sure that in one dimension for s = 1/2 the condition n > 0 is 
sharp for non-existence. 

6.2 Existence for s > 1/2, n < 2s — 1 

Even if the diffusion nonlinearity is singular, the range is formally the same as the “good 
fast diffusion” range 1 > m > (N — 2s)/N, considered in the general theory of [14] , hence 
it is supercritical in the notation of that paper. But we recall that there only exponents 
m > 0 were considered. The presence of the singular nonlinearity when m < 0 makes it 
impossible to use just the same arguments to develop a general existence theory, though a 
number of basic results are expected to be the same. 

We refer to the recent paper |25| by the third author, in which the Barenblatt solutions 
are constructed in this range of parameters; the proofs starts by the existence of a positive 
sub-solution. The sharp behaviour at infinity of Barenblatt solutions is established. 

Finally, we mention that there exist a special class of very singular solutions which are 
solutions of equation ut + (—A ) s (u m /m) = 0 in an appropriate sense specified in [25], and 
has the explicit form 

(6.3) U[x, t) = C (T — l) 1 /^-™) |x|- 2s/(1 “ m) , 

with a constant C = C(m, N,s ) >0 which can be explicitly determined, see [25] . It happens 
that the spatial profile has a non-integrable singularity at x = 0. It is also proved that such 
a formal solution is the limit of a monotone increasing sequence of Barenblatt solutions, we 
refer to [25] for further details. 

The proof of Theorem 12.21 is now complete. 


7 Non-existence for classes of generalized solutions 

The trivial limits obtained in Theorems 12.11 and 12.21 are a clear indication that no solutions 
must exist in any other reasonable sense, since the approximation is monotone from above 
and the maximum principle is known to hold for regular solutions; hence we expect the 
non-existence for other classes of solutions. Indeed, this last fact needs proof, since limit of 
solutions could turn out to be unreliable. 

Loosely speaking, we seek a class of solutions for which the comparison principle holds. It 
turns out that such class is difficult to choose, various non-equivalent classes of solutions 
obey the maximum principle, for example strong or mild solutions. We will prove the 
needed comparison result for strong solutions, because we are interested in short and clear 
arguments. It shall be mentioned that bounded (energy or even very) weak solutions are 
expected to be strong solutions, this is indeed true when m > 0, cf. HUES], but when 
m < 0, this fact requires a different proof, which is quite long and technical and falls out 
the scope of this paper. 
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Theorem 7.1. Let u be a strong solution of Equation (11.11) in the sense of Definition \‘2.1\ 
corresponding to the initial datum uq £ L+(R N ) and let u £ be the solution of the Cauchy 
vroblem \2.3\ . with initial datum u £ $ := uq + e. Then, 

(7.1) u <u £ a.e. in Q = R N x (0, +oo). 


Proof. Let it be a strong solution, hence dtu £ L 1 {Q). We denote with ps a smooth 
approximation of the sign function, namely ps is a a C 1 function with 0 < ps < 1, ps = 0 
for y < 0 and p' s > 0 for y > 0. We take the difference between the equation for u and the 
equation for u £ and we test it with ps(4>[u) — 4>{u £ )). We have 



u - u £ )ps(cf(u) - 4>(u £ ))dx 


[ [p8(<t>(u) - (j>{u £ ))] (-A) s (^(u) 

Jm. N 


cf{u £ ))dx. 


Now, the second integral is non-negative thanks to the Stroock-Varopoulos inequality, that 
reads in this case: 


(7.2) 

[ (ps(4>(u) - <f>{u £ ))) (-A Y((j)(u) - 4>(u £ ))dx > [ 
JR N JR N 


(-a r/^oxu)-^))) 


2 

dx > 0 


where p’ 5 = ( v k^) 2 - We refer to [14], Lemma 5.2 for a proof. Indeed, the above proof is 
formal, since the above inequality holds a class of functions with fast decay at infinity. To 
make it rigorous one should consider further approximations (f{uj)—(j){u Ei j) with a fast decay 
at infinity, for example in the Schwartz class S(R N ) and then pass to the limit, noticing 
that the sign is preserved under such limit process. We leave the details to the interested 
reader. Thus, we have 

(7.3) / dt(u(x,t) — u e (x,t))ps(<j>(u(x,t )) — (f>(u £ (x,t)))dx < 0 for any t > 0. 

J R N 


Since u £ = v £ + e (by definition) and v £ is a strong solution, we have that dtu £ £ L 1 ( l w ) 
for any t > 0. Thus, we can use the dominated convergence Theorem to remove the 6 
approximation in (17.31) . Thus, we get 


d_ 

df 


(■ u(x,t) — u e (x,t))+dx < 0, for any t > 0, 


where we have used that sign + (it — u £ ) = sign + (^>(u) — <f(u £ )) thanks to the monotonicity 
of cj). Thus, integrating the above relation between 0 and some arbitrary t > 0, we get 



u £ (x, t)) + dx < 



dx = 0, 


that implies (17.11) . 


□ 


Clearly, once we have this comparison principle at our disposal the nonexistence result 
applies trivially to strong solutions, thus providing the proof of Theorem 12.31 
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8 Comparison between different nonlinearities 


Up to now we have considered the problems of existence or non-existence of nontrivial limit 
solutions for equations of the type (11.11) in the case of the most typical nonlinearities, those 
of the form f{u) = u m with m > 0, f(u ) = —u~ n for n > 0, of <f>(u) = log(u) (which is a 
kind of case m = n = 0. The results obtained in our previous sections can be extended to 
more general nonlinearities <f using the symmetrization comparison results proved by one 
of the authors and Volzone in papers m, EB- The first one contains the basic comparison 
results for solutions of the same equation, as well as the needed concepts and tools. The 
comparison of solutions of equations with different nonlinearities is presented in m : first, 
we need to recall the concepts of symmetric rearrangement of a positive real function and 
the order relation -< that can be taken from Section 1 of m • Then we need the concept of 
comparison of nonlinearities, also called diffusivities (a bit vaguely). 

Definition 8.1. Assume that (j>, <f> : M+ —> M+ are two functions which are smooth in 
(0, oo). We say T is slower than f if 

$'(r) < <j>'(r), Mr > 0. 

In this case we also say that $ is less diffusive than <f>. 


Once we have introduced the concept of diffusivities , the basic parabolic comparison for 
solutions of the same equation with different diffusivities is proved in 1313 

Theorem 8.1 (Teorem 3.2 [M]). Let u be the nonnegative mild solution to problem 

j u t + (— A) s 4>(u) = / it#, t> 0, 

\ u(x,0) = uo(x) x £ R n , 

with 0 < s < 1, with initial data uq £ L 1 (M Af ), uq > 0, right-hand side f £ L 1 {Q), f > 0, 
and smooth, strictly increasing nonlinearity 4>(u) on with 4>(0) = 0. Assume that is 
a concave nonlinearity on satisfying the same assumptions of <f> and let v be the mild 
solution to 


vt + (—A) s <h(v) = f(\x\,t) x £ R n , t > 0, 
v(x,0) =uq(x) x £ R n , 


where f £ L 1 (Q), uo £ L 1 (M iV ) are nonnegative, radially symmetric decreasing functions 
with respect to x. If moreover d> is slower than 4>, and 


4(\x\) A n 0 (|x|), f*(\x\,t) ~< f(\x\,t), 


for almost all t> 0, then the conclusion u#(\x\,t) v(\x\,t) holds. 


We apply the Theorem above to a pair of nonlinearities <f(u) and 3> n (rt), where $> n (u) = 
—u~ n for n > 0 and 3 >q(^) = log (it). We choose the case where n satisfies the conditions of 
the non-existence theorems, and we get the following improved result. 
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Theorem 8.2. The nonexistence results of Theorem 12.11 and Theorem m still hold when 
instead of taking <f> = we assume that cj) is smooth, strictly increasing nonlinearity <f{u) 
on M_|_ with (f>( 0) = 0, and moreover 

(8.2) <j)'{u) > C(M ) u~^ +1 ^ for all 0 <u< M, 

and all M > 0. The range of accepted ( N,s,n ) is as in Theorem \2.1\ and Theorem \2.2[ 

Corollary 8.3. If a non-existence result is proved for a choice of (. N , s, n), it also hold for 
(N, s, n i) with n± > n. 

This corollary sheds light into the conditions of the phenomenon of non-existence. 


9 The elliptic problem 

The techniques used to study the parabolic problem allow us to prove analogous non¬ 
existence results for the following elliptic equation 

(9.1) u + (-A) s cj){u) = f in R n , 

when / € L 1 (M iV ) and the nonlinearity 4> is given by the formulas (12.21) . Writing f>(u) = v 
hence u = (f~ l {v ) := /3(v), we can rewrite it in the form 

(9.2) (-A ) s v + /3(v) = f in R N . 

This is the fractional version of the equation treated by Benilan, Brezis and Crandall, for 
s = 1 in their famous paper [2]. In the paper, they study existence for equation (19.21) 
assuming that /3 is a monotone function. More precisely, they assume that /3 is a maximal 
monotone graph (hence possibly multivalued) with the property that 0 € /?(0). On the 
other hand, our nonlinearity satisfies j3(— oo) = 0, which is the cause for the non-existence 
result that is not considered by them. 

Equation (19.11) appears as the iteration step on solving the parabolic equation (11.11) by 
the implicit time discretization scheme suggested by the Crandall-Ligget Theorem [9j. This 
explains why the solutions we get for the elliptic problem, Theorems 19.11 and 19.21 are so 
similar to the parabolic ones. 


9.1 Nonexistence of limit solutions 


We consider an approximate problems with / G L 1 (M jV ) D L°°(Ml n ) 


v + (-A ) s ((j) £ (v)) = f in 1 A ' 
4> £ {v) := 4>(v + e) -4>{s). 


Notice that for any e > 0 the nonlinearity cj) e is smooth and monotone increasing. Therefore, 
the above problem can be solved in a standard way (see the discussion in j24l Section 8.1]), 
and the solution turns out to be unique, bounded and nonnegative. Moreover, 


(9.4) 


\\ v c\\l° 


< 111 llL° 
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and the standard comparison theorem holds. Let us denote by v £ such solutions. Next we 
define u £ := v £ + e, noticing that u £ is the (unique) solution to 

(9.5) u £ + (-A) s 4>(u £ ) = f £ := f + e in . 

We also have 

(9.6) u £ >e in R N , 

(9-7) ll^e ^■IIl 1 (r jv ) — II/IIl 1 (e jv ) 

Recall that the sequence u £ is monotone with respect to e, therefore the limit function 

(9.8) u(x) := limu e (x) 

exists for almost any x G M. N . Moreover, we have that 0 < u(x) < ||/||l°o( K jv). In the case 
where / G L 1 (M iV ) we can use the property of L 1 contraction to show that there is also a 
limit for formula (19.811 though the function u £ are not necessarily bounded. 

Analogously to the parabolic case, we will prove the following result. 

Theorem 9.1. For any given f with f G L 1 (M jV ), let u £ be the solution of (19.51) . Then the 
trivial limit 

(9.9) lim u £ = 0 in ^(l^) 

e\0 

holds true whenever N > 2 and n > 0, or N = 1 and 

n > 0 and s G (0,1/2); or n > 2s — 1 and s G (1/2,1); or n > 0 and s = 1/2. 

This range of exponents is optimal. 


Proof. The proof will be split into several steps. We repeat the scheme of the proof for 
the parabolic case, therefore we will just sketch the main arguments. We will begin with 
/ G L 1 (M JV ) D L°°(R n ) and in Step 5 a density argument will allow us to treat the case 
/ G L 1 (M Ar ). The last step is devoted to prove the sharpness of the exponent’s range. 

• Step 1. The case N > 2s. This case includes the proof of both the case N > 2 and 
N = 1, s G (0,1/2). Proceeding as in Lemma 14.21 multiplying by the test function 


V'n(') 


(-A)- s 


XB 1/n (x) \ 


and passing to the limit n —> oo we obtain 

(9.!°) <Ku e (x)) - 4>{e) = J^ n £ dy. 


Notice that the right-hand side is uniformly bounded with respect to e > 0 and x G M. N , 
because 

(9-11) [ /( f - U |ff 2a +£ dy < 2Cmax{||/|| Ll(RJ v ) , ||/|| i «» (RJ sr ) } 

jr n f — y\ 
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the proof being identical to the parabolic case, see formula (14.141) in the proof of Lemma 
IQ As a consequence, the following limits are uniform in x € 


lim 4>(u £ (x)) = —oo , lim u e (x) = 0 . 

e\0 e\0 


This clearly implies that u = 0 in 

• Step 2. This case corresponds to N = 1, s € (1/2,1), and n > 2s — 1. Let 

u(x) := lim u £ (x), 
e\0 

and assume by contradiction that u( 0) > 0. Proceeding as in Lemma 15.11 we obtain the 
inequality 



(9.12) 


Thus, in the limit £ \ 0, the same inequality holds for the limit function u that is 



Finally, we obtain a contradiction since n > 2s — 1 and u € L 1 (M) by construction. 

• Step 3. This case corresponds to N = 1, s = 1/2 and n > 0. We use the same strategy 
of the corresponding parabolic problem. For any e > 0 we approximate the solution of 


u £ + (-A) 1 ' /2 ^(u e ) = f £ 
fs = f + e 


(9.13) 


with solutions of 



problem, cf. m Theorem 10]. Hence, for S small enough, we have that u £j s satisfies the 
lower bound (19.121) . namely 



Next, we let <5 \ 0 and recall that the constant C does not depend on s neither on 5, hence 
we have that u £ satisfies the desired lower bound 
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As a result, we have a contradiction once we assume that the limit function u = lim £ \ j o u e 
is not identically equal to zero. 

• Step 5. To conclude the proof, we need to discuss the general case of unbounded forcing 
functions in L 1 (R Ar ). This can be done exactly as in the parabolic case, through a similar 
density argument. 

• Step 6. Sharpness of the range of exponents. We only have to show that the limit maybe 
nontrivial for some exponents s and n is one space dimension. Using the idea of paper [25] 
for the parabolic problem, we only need to show that there exists some nontrivial solution 
or subsolution that lies below the approximations u e . 

• In the borderline case s = 1/2 and n = 0, we take the spatial profile of the explicit solution 

m, 


F(X) = y 


+ \x\ 


that satisfies (-A) s (log(i ? )) = 2 F in R. Therefore F is a nontrivial and integrable solution 
of the equation 

u + (-A) s (log(ri)) = f =3F . 

• In the case s > 1/2 and n > 0, n < 2s — 1, we use a similar approach using a special 
subsolution constructed in Section 3.2 of ESI. □ 


9.2 Comparison between solutions and approximate solutions 

We want now to prove non-existence for a class of weak solutions, defined as follows. 
Definition 9.1. A function u is a weak solution of equation (19.11) if it is nonnegative and 

• u G L 1 (R N ), <t>(u) G H s (R n ) 

• It satisfies, for all £ G Cq°(M' v ), 

(9.15) f u(x)if(x)dx + [ (-A) 5 / 2 (j>{u{x)) (-A) S//2 if(x)dx = f f(x)'if(x)dx. 

Jr n Jr n Jr n 

As in the parabolic problem, our goal now is to prove the following non-existence result 
for the above class of weak solution to the elliptic equation m, which is the elliptic 
counterpart of Theorem 12.31 

Theorem 9.2. Let u denote a nonnegative weak solution of the elliptic equation m, 
corresponding to some f G in the range of parameters ( N,s,n ) allowed in Theorem 

19.11 Then, u = 0 in W N , hence f = 0. This range of exponents is optimal. 


The proof of the above Theorem follows by combining the results of Theorem 19.11 for 
approximate solutions together with the following comparison result between approximate 
solutions and weak solutions, which is the elliptic counterpart of Theorem 17.11 
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Proposition 9.3. Let f be a given L 1 (M jV ) function and let u be a weak solution to the 
elliptic equation m in the sense of Definition 19.11 Then, 

(9.16) u < u e a.e. inW N , 

where u £ is the solution (19.51) with right hand side given by f + e. 


Proof. The proof is similar to the proof of Theorem 17.11 therefore we just sketch it. We 
take the difference between the equation solved by u and the equation solved by u £ and 
then we test the equation with pfifi(u) — fi{u £ )), with ps the same smooth approximation 
of the Heaviside function used in the proof of Theorem 17.11 We have 



u - u £ ){x)p s (<j)(u) - f(u £ ))(x) dx 


+ 



[ps(<t>{u) - 4>(u £ ))] (-A) s (<f(u) - 
(-e)ps((j>(u) - c/>(u £ ))(x)dx < 0. 


fi(u £ ))dx 


The second integral in the left hand side is positive thanks to the Strook Varopoulos in¬ 
equality CL2D- On the other hand, since u and u £ are in L 1 (R JV ), we can let 6 \ 0 and 
obtain 

[ (u- u £ )+{x)dx < 0, i.e. u < u £ for a.a. x € 

Jr n 

where we have used that sign + (u — u e ) = sign + (0(u) — cj>(u e )). □ 


10 Non Existence for the Dirichlet problem 

In this last Section we discuss the question of existence/non existence for solutions of a zero 
Dirichlet problem for (12.11) in a bounded domain O C WL N . In this case, the scenario for 
existence is identical or worse than the Cauchy problem in the whole space. In particular, 
for standard diffusion (i.e. s = 1) solutions do not exist for integrable initial conditions in 
any dimension if m < 0 (see [28] 1. 

We recall that the fractional Laplacian can have different nonequivalent expressions on 
bounded domains, hence we have to study two different Dirichlet problems. There are two 
choices of fractional laplacian on a domain are: the Spectral and Restricted Laplacian , cf. 
for the definitions. Then, the two different Dirichlet problems are as follows 


(10.1) 

f dtu + (-A ) s ((f>(u)) = 0 in P x (0, +oo) 
[u = 0 mR N \Q, 

or 


(10.2) 

[dtu + C s {(j){u )) = 0 in 9 x (0, +oo) 

= 0 on dD. 
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where we have denoted with C s the Spectral laplacian (see 0 )- 


Independently of the choice of the type of fractional laplacian, we have the following non¬ 
existence result 

Theorem 10.1. Let Li be a bounded domain in R N . Let u denote a nonnegative strong 
solution of (110.111 or of (110.21) in the range of parameters ( N,s,n ) allowed in Theorems \2.l\ 
and a Then, u = 0 in R N x (r, +oo) for any r > 0. 


The above Theorem will be proved using a comparison principle. In particular, it is 
known that for smooth nonlinearities (like our </> e ) the solutions of the Dirichlet problem 
for the Spectral laplacian lie below the solutions of the Cauchy problem, cf. e.g. mm, 
hence by comparison we deduce the same non existence results of Theorem 12.31 also for the 
homogeneous Dirichlet problem. 

Here we discuss the non existence for the problem with the Restricted Laplacian. As 
expected, the proof of the Theorem is based on a comparison argument for the approximate 
solutions of (12.111 and (IIP.111 . More precisely, consider rto £ L 1 (R N ) and fix some compact 
domain Q' C D. Then, for any e > 0 we denote with w e the solution of 


(10.3) 


d t w £ + (-A) s cj) £ (w E ) = 0 infix (0, +oo), 
< rc e (x,0) = xn'{x)u 0 (x) in R N . 
w E = 0 in (R N \ Q) x (0, +oo), 


where, as in the previous section, 4> £ {y) := f{y + s) — 4>(s)- We recall that the Restricted 
fractional Laplacian is defined, for a smooth function v supported in D as 


(-A) s v{x) := c(N , s)p.v. f iy, for x G R N . 

j«. n \ x — y i 

The well posedness for (110.31) follow from flljllSllSKllE]- Incidentally, note that w E := w £ +e 
solves 

dtw £ + (-A) s 4>(w e ) = 0 in fi x (0, +oo), 

< w £ (x, 0) = xn'(x)uo(x) + e in R N , 

^ w £ = £ in (R n \ D) x (0,+oo). 

Moreover, for all points in (R N \Q) x (0,oo), we have that w £ (x,t ) = f> £ (w £ ) = 0, hence 


d t w £ (x,t ) + (-A ) s (f £ (w £ (x,t)) = c(N,s) p.v. 

Jr 1 


(j) £ (w £ (x,t)) - <f> £ (w E (y,t)) 


\x - y | 


N+2s 


d y < o, 


where cj) £ is a positive function. Consequently, we have that w £ verifies 


dtw £ + (-A) s (j) £ (w £ ) <0 in x (0, +oo), 


with equality when x G D. Thus, since Xn lU o A uo, the comparison principle gives that 
We ^ v e in x (0, +oo) (hence w £ < u £ ) where v £ is the solution of (13.11) with uq as initial 
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condition. On the other hand, it is not difficult to show that any weak solution of 


dtw + (-A) s <j>{w) =0 in 0 x (0, +oo) 

< w = 0 in (R N \ D) x (0, +oo) 
w(x, 0) = xw{x)uq{x) in R N , 

satishes w < w £ in Q. Thus, we have the chain of inequalities w < w £ < u £ and hence by 
letting e \ 0 and recalling Theorems 12.11 and 12.21 we conclude. 


11 Comments and open problems 

• Our result for the Dirichlet problem may not be optimal in dimension one. In fact, since 
our proof (both for the Spectral and for the Restricted Laplacian) follows by a comparison 
argument, we actually prove that the non-existence range for the Dirichlet problem is at 
least as large as the non existence range for the Cauchy problem. Is it larger? this question 
remains an open problem. There is some evidence that may help. On the one hand, in 
the cases (N = 1 ,s = 1/2, n = 0) and (N = l,s > 1/2, n < 2s — 1) the Cauchy problem 
has a solution (see Section [6] in this paper and the recent [25]). On the other hand, for the 
Dirichlet problem in the case N = 1 and s = 1 (namely, standard Laplacian) non-existence 
holds for any exponent n > 0. 
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